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• In this talk we consider a non standard evolution for the B'^-B^system, namely, an evolution in 

On ! the open quantum systems framework. Such approach is justified by the very high sensitivity 

' of experiments studying CP-violating phenomena in the B-mesons sector, very near to the one 

required to test some possible scenarios induced by quantum gravity at the Planck scale, whose 
effects at low energy can be described by a heat bath. We adopt a phenomenological point of 
view, introducing six new parameters that fully describe this kind of evolution without referring 
to a specific model for the microscopic interaction. We outline the main differences between 
Cl^l this approach and the usual one in the description of evolution and decay of single mesons or 

^ • correlated pairs. 



1 Introduction 

The physics of b-flavoured hadrons is at present one of the most promising fields testing Standard 
Model predictions concerning some interesting physical processes like CP-violating phenomena 
and flavour oscillations. Presently, several experiments are collecting data; other very accurate 
experiments are planned to start in the near future. In this context it is also possible to study some 
exotic scenarios of physics beyond the Standard Model; in this contribution we consider an open 
quantum system treatment for the 6*^-6 system P producing decoherence^ . 

A wide variety of physical systems has been studied in this framework. The formalism of open 
quantum systems, initially developed in quantum optics |3ISlin|, has been also applied to study 
phenomena of quantum relaxation in magnetic resonance [3, to foundational aspects of quantum 
mechanics [5] and, lately, to model some typical situations in particle physics: neutral mesons 
systems (HI , neutrinos J2] , photons ^Hl , neutrons interferometry jll] . 

An open system |151 1161 117j is a system S that interacts with its external environment E, and, 
then, exchanging energy and entropy with it. Whereas the dynamics of the composite system 
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^For a preliminary analysis see 0; for a different approach to decoherence in the B-mesons system see 
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T = S + E is unitary, this is no longer true for the dynamics of the system S alone, obtained 
tracing the time-evolution of T over the degrees of freedom of the environment E. In general, the 
resulting reduced dynamics is complicated by memory effects. It is possible to get a Markovian 
evolution for the system S (that is, a memoryless dynamics) if essentially two assumptions hold: 
the interaction between S and E must be weak, and there must be no initial correlation between 

them [miiniini. 

Since we are interested in time-evolutions causing decoherence, that is transitions from pure 
states to mixtures, we must describe the state of S by means of statistical operators p, that is 
Hermitian and positive trace-class matrices with unit trace. The time-evolution for p is given by 



p{t)=LW)] = -iHp{t) + ip{t)H^ + LdW)], (1) 

where L is the generator of the dynamics, H is the Weisskopf-Wigner Hamiltonian describing the 
system S (non-Hermitian because the system is unstable) and is the generator of the dissipative 
part, emboding the interaction with the environment. The form of Lo, univoquely fixed by the 
Kossakowski theorem |15| I16j.is very general since it characterizes any Markovian dynamics: 

LdW)] = -\[Y,A\A,,p{t)}+Y,A.^p{t)A\, (2) 

i i 

where the set of operators Ai depends on the details of the interaction. 

The time-evolution maps /3(0) — > p{t) = At[/9(0)], where At = e^*, satisfy the forward in time 
composition law: A^+s = A^ o A^, with t, s > 0, characteristic of irreversible phenomena; further- 
more, they are completely positive |15l I16j . This last property is noteworthy in particle physics, 
because it guarantees the preservation of positivity of evolving density matrices describing entan- 
gled systems. Indeed, simple positivity is not sufficient for a proper statistical interpretation of the 
formalism of density matrices when they represent entangled systems |151 116j (counter examples 
are found in ^2)- Moreover, the evolution maps At allow transitions from pure to mixed states. 

The entropy of open systems can increase as well as decrease; nevertheless, in what follows, we 
will consider only entropy increasing time-evolutions because we are interested in reduced dynamics 
producing states less ordered as time passes. 

As we shall shortly see in the next section, the phenomenonogical approach to the B'^-B^system 
is usually performed in the two dimensional Hilbert space of flavours; then we shall consider 2x2 
density matrices p that can be expanded over the Pauli matrices fJi [i = 1,2,3) and the 2x2 
identity ctq: 

3 

P = ^Pi(^i- (3) 

i=0 

It is possible to give a vectorial representation of the density matrix p, namely, defining a 4- vector 
\p) whose components are the real coefficients pi in Q: 



Pi 

P2 
\P3 J 
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accordingly, since equation is linear in time, it can be rewritten as a Schrodinger-like equation: 



dMt)) = {n+v)\p{t)), 



(4) 



where Ti. and T> are 4x4 real matrices, the Hamiltonian and the dissipative part respectively. The 
non-standard part of the dynamics is embodied in the symmetric matrix D: 



V 



/O 0\ 

a b c 

Q h a (3 

VO c /3 7/ 



(5) 



with a, 6, c, a, (3, 7 six real phenomenological constants, constrained by the inequalities, necessary 
and sufficient to Complete Positivity [T3] : 



2R = a + -f - a > 0, RS>b^, 
25 = a + 7-Q > 0,RT>c^, 
2r = a + a-7 > 0, ST>(3^, 
EST > 2bcp + + 5c2 + r62. 

Since these coefficients are smalP with respect to the entries of Ti., it is possible to solve 
to first order in these parameters: 



(6) 



up 



\p(t)) = e^'m) + fdse^^'-^^Ve^^pm- (7) 
Jo 

The first term describes the standard evolution of the density matrix, the second one accounts 
for the interaction with the environment. 

In the following we won't report the explicit form of statistical operators at time t because it 
is not physically relevant. The computations have been performed in vectorial representation and 
using the approximation Q. 



2 Decay of neutral B-mesons 

The neutral B-mesons system can be represented by statistical operators p acting on the two 
dimensional flavour Hilbert space, spanned by the physical states IB'') and |B ). In a dissipative 
context, their time-evolution in vectorial representation is given by ((7j); the computation has been 
performed working in the basis of the eigenstates of the Hamiltonian H, denoted by |Bl) and |Bh), 
where the subscripts mean "light" and "heavy" respectively: 

i7|BL) = Az.|BH), ' ^^'L-mH,L-^lH,L (8) 



rough evaluation of their magnitude can be done: Vij ~ e|/e_B, where eg and es are the characteristic energies 
of the system and of the environment respectively 
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where rriH^L and '~ih,l are the masses, respectively the widths, of the corresponding eigenstates. 

This choice of basis provides an invariant formahsm under independent phase redefinitions of 
the two basis vectors. Therefore it makes no commitment on the magnitude of CP (and CPT) 
violating effects. 

In the density matrix formalism, the mean value of an observable O is obtained from: 

{0{t))=Tr[Op{t)]. (9) 

Then, the probability rates for the decays of B'^and B'^in different channels, B", b''— > /, can be 
written as 

Vf{B^,t) = Tipfpit)], with p(0) = |B0)(B0|; 
Vf{^,t) = Tv[Ofp{t)l with p(0) = |B°)(B°|, 



where the operator O/ is defined by: 



with entries: 



Of=r (10) 



01 = \A{^^ f)\\ 03 = [A{^^ /)]M(b"^ /), 

02 = |^(B°^ 04 = [^(B°^ /)]M(b0^ /), 

where A the the decay amplitude for the considered process. 

The standard dynamics is characterized by Am = mn — A F = 7^ — 7// and 

r = 7i + 7_f//2; it is customary to renormalize these parameters, and the time variable, in the 
following way: 

Am AT , , 

^ = — , ^ = ^^ ■^ = *r. (11) 

The dissipative contribution is instead parametrized by particular functions of a, 6, c, a, (5 and 7, 
denoted by A, C and D; their expressions are reported in 



2.1 Decay of a single B-meson 

Typical decays, relevant for the study of CP (and CPT) violation, are the so-called semileptonic 
decays, in which the final state / contains a hadron h, a lepton / and a neutrino (or anti-neutrino): 
B0(B°)^ h-l+v{h+l-V). 

It is customary to parametrize the amplitudes in (fTU]) as follows ^1 : 

^(BO^ h-l+u) = Mh{l-yh), 

A{^^h+l~u) = Ml{l + yl), 
^(BO^ h+l^v) = ZhA0^^ h+l^v), 
^(B°-^ h~l+u) = XhAiB^^ h"l+u), 
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where the parameters Xh and Zh measure violations of the rule AB = AQ, satisfied in the Standard 
Model, and the coefficient yh measures deviations from CPT invariance. 

For sake of simplicity, the results we show below are obtained assuming Xh = Zh = yh = 0; more 
complete results are collected in The probability rates present new dissipative terms absent in 
the standard case; for example, in the case of the decay B*^— > h^l~V: 



Vh+{B\T) = ^^^e-"|cosh5r-coswT e-(^"^)"+ 
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+ sinh5rf^ - -j^^Im(C)') + sincjrf-r^^ Im (C) - Re 



(12) 

If dissipation is neglected (A = S = C = = 0) we obtain the standard expression of the 
probability rate: the second row of H12|l vanishes and the damping exponential in the first row 
disappears. 

From a fit of the time dependence of (|12|) it should be possible, at least in principle, to get the 
parameters yl, Z?, Re(i?) and Im(C). 

Usually, instead of analyzing directly the decay rates, some particular asymmetries, constructed 
by means of these rates, are taken into account. The following one, for example, is extensively 
studied: 

Indeed, in the standard case, this quantity is used to fit the parameter Am introduced previously; 
when we consider the non standard evolution it assumes the simple form: 

-4Am(T) = e~'^'^cosa;T + Re (i?) sincjT (14) 

that reduces, in the standard case (i.e. A = B = C = D = Q) to ^Am('7") = coswr. We see 
that a study of the time dependence of this quantity should provede, in principle, the dissipative 
parameters A and Re(-B). 

The evaluation of the parameter Am in the standard case is also performed considering time 
integrated rates 

Vf{B) = ^ dTVf{B,T) (15) 

and introducing the asymmetry defined in analogy to (|1.S|) . with the time integrated rates; 

the expression of this quantity in the dissipative context is: 

-4'a™ = yT:J2 { 1 + 5? (^) + [(^' - 1) ^ - 2a;'^] } • (16) 

In this case, since the dependence on time is lost, it is more difficult to constrain the parameters 
A, Re(i?) and D. This is a general problem for time-integrated quantities. 
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Another class of relevant decays of single B-mesons has a final state that is a CP eigenstate, for 
example: / = D"^D~, 7r"^7r~, i/ip K. In this case it is convenient to define a new asymmetry: 

This quantity would be rigorously zero in absence of CP violation; it is studied because it enables 
to fit the parameter P appearing in the CKM matrix |18l I19j . With the non-standard contributions 
to the time-evolution it takes the form: 

Afit) =^Im (C) - sinu^r ( ^^^^ + ^Re {C)) + 

U + IA/P J ^^^^ 



cos LOT ( ^ - ^Im (C) 

1 + Af r uj 



where A/ is a phase- independent coefficient introduced in order to parameterize the neutral mesons 
decays ^lE]- A fit of this asymmetry would provide information about the parameter C. 

2.2 Decay of correlated B-mesons 

Of particular interest are the decays of pairs of entangled B-mesons, produced at the so-called 
B-factories. As already stated, the property of Complete Positivity plays a fundamental role in this 
case. 

The entangled B-mesons are generated from the decay of the T(4S) resonance: T(4S)^B'^B ; 
the initial state is Ii/jb) = (I ^'^i ~P) ® — jB*^, —p) (g) jB", p)), and the initial density matrix 

p{0) = \iIjb){'4'b\- Its time-evolution is supposed to be factorizable: 

p(ti,t2) = (At, ®AjJp(0), (19) 

where the first meson has evolved up to a time ti and the second up to t2. The probability of decay 
for the first meson in a final state fi at time ti and for the second meson in /2 at t2 is given by 

V{h,h;f2,t2)=Tr[{Of,®Of,)p{h,t2)], (20) 

where the operators Of have been introduced in (|1U() . 

In the experimental analysis single time distributions are usually taken into account; writing 
ti = t' + t and t2 = t' we express these probabilities by 

r+oo 

^/i,/2(0= / dt'V{fi,t' + t;f2,t') (21) 
Jo 

for positive t, i.e. when the first meson decays after the second one, and by 

r+oo 

'Pfuf2i-\t\) = / dt'Vih,t' + t;f2,t')eit'-\t\) (22) 
Jo 
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for negative t (when the first meson decays before). These definitions fulfill the logical request: 

'Pfuf2i-\t\)=rf,jM)- 

We have analyzed the behavior of some characteristic quantities of the double semileptonic 
decays, i.e. those decays in which both mesons undergo a semileptonic decay.The first quantity we 
consider is the asymmetry: 



Its relevant property is that in the standard case 7^(0) = 0, because 'P/,/(0) = 0, whereas taking 
in account the dissipative contributions it becomes 



^(0) = ^ 



A + 



1+UJ' 



(wlm(5) - Re (B)) 



(24) 



regardless of the assumptions made about the coefficients Xh, Zh and yh- Then a non-zero value of 
7^(0) clearly shows the presence of dissipative phenomena. 
Another interesting asymmetry is expressed by 



. . . [Ph+,h-{r)+rh'Mi^)] - [n+,fe+(T)+n-.fe~(T)] 

It is analogous to the one defined in for the single meson system, because it allows, if there 
is not dissipation, the determination of the parameter Am. In our case the dependence on the 
dissipative parameters is: 

e~^'^ B 
AAm(T) = 7 coswr + (a; coswr + sinwr) Re ( ). (26) 

1 + A 1 — ILO 

In the above expression we have assumed 6 = 0, a good approximation since AT <C F.From 
the time-dependence of the contributions in H26|) it should be possible to evaluate the dissipative 
parameters A and B. 



3 Conclusions 

We have treated the neutral B mesons system in the framework of quantum dynamical semigroups. 
The formalism is very general and can be applied to describe the time evolution of this system when 
it is subject to a weak interaction with an external environment. For instance, quantum gravity 
could give a motivation for such an approach, since in this case the space-time itself would act as 
an effective environment for any physical system. 

The open system framework allows a rough evaluation of magnitude of the dissipative parameters; 
they scale at most as m^/Mp, where me is the meson mass while Mp is a large fundamental mass 
scale. Tipically, Mp coincides with the Planck scale so that m^/Mp ~ lO"'^^ GeV. This value 
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is very small; however, the sofistication of the dedicated B-experiments, both at colliders(CDF-II, 
HERA-B, BTeV, LHC-6) and B-factories (BaBar, Belle, CLEO-III), is so high that the sensitivity 
required to probe the presence of non-vanishing dissipative parameters can be reached in just a 
few years of data taking |2n | I2 H [22]. If these parameters should be found not different from zero, 
the inequalities expressing the Complete Positivity could be tested. Then, the neutral B mesons 
physics is potentially a good laboratory to test basic properties of open systems dynamics. 
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